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ABSTRACT 


■  The  resolution  of  an  in-line  hologram  recorded  with  plane  waves 
in  the  Fraunhofer  region  of  a  circular,  opaque,  particle  is  evaluated. 
From  the  diffraction  integral  for  the  reconstructed,  real  image  the 
intensity  distribution  about  the  image  is  found.  Criteria  for  deter¬ 
mining  the  edge  of  the  image  in  addition  to  determining  the  plane  of 
focus  are  specified  from  the  data,  and  the  measurement  accuracies 
are  found.  The  inaccuracies  are  shown  to  be  a  result  of  the  inabil¬ 
ity  of  the  film  to  record  all  of  the  light  diffracted  by  the  particle. 

The  reasons  for  this,  film  grain  noise  and  its  cutoff  frequency  and 
dynamic  range,  are  explained  and  their  relative  effects  compared  for 
the  various  particle  sizes  and  recording  distances. 
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1.1  PROBLEM  STATEMENT 

In-line,  Fraunhofer  holography  is  a  convenient  means  of  recording 
and  measuring  particle  size,  position,  and  velocity  in  a  dynamic  volume. 
Over  the  past  few  years  it  has  been  used  in  aerosol  studies  to  measure 
the  size  distribution  and  velocity  of  the  particle  droplets  (Refs.  1,  2, 
and  3).  These  studies  have  brought  out  the  need  to  know  the  limiting 
accuracy  for  measuring  the  size  and  position  of  the  reconstructed  images. 
The  accuracy  of  these  measurements  is  important  especially  for  the 
velocity  data  where  the  ability  to  bring  an  image  into  focus  and  determine 
its  exact  position  from  the  hologram  is  essential.  The  factors  which 
limit  this  accuracy  must  also  be  known  so  that  the  hologram  recording 
may  be  improved  to  yield  more  accurately  reconstructed  images.  The 
errors  inherent  in  measuring  the  particle  image  size  and  finding  the 
planes  in  which  the  reconstructed  images  are  in  focus  should,  therefore, 
be  determined  analytically. 


1.2  IN-LINE  HOLOGRAPHY  AND  ITS  APPLICATION  TO  INSTRUMENTATION 

There  are  two  basic  ways  of  recording  a  hologram:  (1)  in  line  and 
(2)  off  axis.  The  simplest,  and  historically  the  first,  is  the  in-line  holo¬ 
gram.  It  is  formed  by  illuminating  a  volume  consisting  of  particles  with 
a  coherent  beam  of  light.  The  light  diffracted  by  the  particles  interferes 
with  the  undiffracted  portion  producing  an  intensity  variation  which  contains 
both  the  phase  and  amplitude  information  of  the  light.  These  intensity 
variations  are  recorded  as  interference  patterns  on  film.  The  developed 
negative  is  a  hologram  which,  when  placed  into  a  beam  of  coherent  light, 
will  reconstruct  a  real  and  virtual  image  of  the  particle  field.  If  plane 
waves  are  used  to  record  and  reconstruct  the  hologram,  the  images  are 
reconstructed  exactly  to  size  and  position  as  the  original  particles  were 
relative  to  the  film.  The  images  are  located  on  either  side  of  the  film 
"in  line"  with  the  reconstruction  beam.  The  recording  and  reconstruc¬ 
tion  arrangements  are  illustrated  in  Fig.  1,  Appendix  I.  The  hologram 
is  recorded  in  the  far-field,  Fraunhofer  region  of  the  particles  to  keep 
the  two  images  from  interfering  during  the  reconstruction  (Ref.  4). 

Sizes  and  positions  of  the  particles  in  the  volume  can  be  determined 
from  the  hologram.  The  velocity  of  the  particles  can  also  be  found  in  two 
ways.  If,  during  the  exposure  of  the  film,  the  particles  move,  their 
reconstructed  images  will  be  smeared.  (Unlike  the  other  holographic 
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techniques,  motion  during  the  recording  is  not  limited  to  less  than  a 
wavelength. )  From  the  length  of  this  smear  and  the  exposure  time  the 
velocity  of  the  particle  can  be  calculated  (Ref.  5).  The  second  method 
of  obtaining  the  velocity  is  based  on  doubly  exposing  the  hologram.  The 
holograms  taken  at  two  instants  of  time  reconstruct  volumes  containing 
similar  particles  that  are  displaced  in  space  according  to  their  velocities. 
Velocities  are,  therefore,  found  from  the  separation  of  the  two  particle 
images  and  the  time  separation  between  exposures  (Refs.  3,  6,  7,  and  8). 

Besides  reconstructing  the  images  and  measuring  their  sizes, 
positions,  and  velocities  in  the  volume,  there  are  two  other  methods  of 
reducing  these  data.  The  first,  and  most  basic,  is  by  direct  analysis 
of  the  recorded  interference  patterns  (Refs.  7  and  9).  This  was  used  in 
the  preliminary  stages  of  particle  field  analysis  using  coherent  illumina¬ 
tion  and  has  been  found  to  contain  a  number  of  shortcomings,  making  data 
reduction  impractical  (Ref.  10).  Optical  data  processing  of  the  hologram, 
the  second  method,  is  as  yet  in  its  preliminary  stages  of  development 
(Ref.  11). 

To  produce  a  high  quality  in-line  hologram,  the  recording  must  be 
made  in  the  far-field  region  of  the  individual  particles.  This  condition 
limits  the  maximum  particle  size  that  can  be  recorded  since  the 
Fraunhofer  region  of  particles  larger  than  one  millimeter  is  more  than 
one  meter  away.  Diffraction  patterns  of  particles  less  than  10  pm  are 
difficult  to  record  because  of  the  small  amount  of  light  diffracted.  The 
density  of  the  particles  in  the  volume  must  not  be  so  large  that  the 
reference  beam  is  cut  out.  Reference  10  gives  a  plot  of  the  maximum 
allowable  density  for  a  given  particle  size  and  volume  depth. 

The  velocity  that  can  be  recorded  is  limited  by  the  pulse  width  and 
separation.  For  nonsmeared  reconstructions  in  the  double-pulse  tech¬ 
nique,  the  pulse  width  must  be  such  that  the  particles  move  less  than  a 
tenth  of  their  diameter.  The  pulse  separation,  however,  must  not  be 
so  close  that  the  images  cannot  be  resolved  or  so  far  apart  that  the  two 
images  cannot  be  paired. 

The  accuracy  with  which  the  particle  size  and  velocity  data  from  the 
reconstructed  images  of  the  particles  can  be  measured  must  be  known  in 
order  to  determine  the  inherent  error  of  the  reconstruction  process.  The 
theoretical  error  involved  in  an  instrumentation  readout  of  these  image 
parameters  limits  the  accuracy  of  the  instrument.  The  study  to  be 
reported  here  was  concerned  with  the  hologram  resolution  of  the  particle 
image  and  its  focal  plane.  Associated  with  every  hologram  is  a  limiting 
aperture  which  affects  its  resolution  capabilities.  The  causes  of  the 
limiting  aperture  and  the  resolution  capabilities  of  a  given  film  are  shown. 
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Section  II  will  consist  of  the  derivation  of  the  focusing  integral  for  the 
image  and  special  cases  of  this  integral  which  have  closed  form  solu¬ 
tions.  In  Section  III  the  data  from  the  numerical  integration  of  the 
focusing  integral  will  be  used  to  define  criteria  to  describe  the  capa¬ 
bility  of  an  instrument  in  determining  the  particle  image  size  and 
location.  The  limiting  criteria  are  found  in  terms  of  the  finite  size 
of  the  recorded  interference  patterns.  The  importance  of  the  various 
causes  of  the  limited  size  of  the  patterns  are  studied. 


SECTION  II 

THEORETICAL  ANALYSIS 


2.1  DERIVATION  OF  THE  FOCUSING  INTEGRAL 

An  in-line  hologram  can  be  recorded  with  plane  or  spherical  waves. 
To  keep  the  analysis  from  being  restricted  to  any  special  case,  the  initial 
derivation  for  the  reconstructed  image  will  be  for  ispherical  waves.  The 
point  source  and  recording  plane  will  be  in  the  far-field  region  of  the 
circular  particle  where  the  recording  distance,  zj,  written  in  terms  of 
far  fields  is 


where  is  the  wavelength  of  the  illuminating  wavefront, 'N  is  the  far- 
field  number,  and  a  is  the  particle  radius.  When  this  criterion  is  not 
met,  the  film  is  in  the  Fresnel  region  and  the  reconstructed  image 
suffers  from  interference  with  the  still  large  intensity  conjugate  image. 
The  schematic  arrangement  for  recording  the  diffraction  patterns  and 
reconstructing  the  images  is  shown  in  Fig.  1. 


The  normalized  intensity  distribution  in  the  recording  plane  caused 
by  an  opaque,  circular  particle  illuminated  with  ^spherical  waves  of  radius 
zQ  is  (Ref.  12) 


I(r) 


(2) 


This  intensity  distribution  is  the  interference  between  the  diverging 
background"light  which  has  been  unaffected  by  the  particle  and  the  light 
scattered  by  the  particle.  The  first  two  terms  represent  a  background 
intensity  caused  by  the  reference  and.  object  waves  and  give  no  informa¬ 
tion  to  the  reconstructed  image.  These  terms  merely  represent  the 
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intensity  that  would  be  recorded  if  the  particle  were  illuminated  with 
incoherent  light,  i.  e,  ,  interference  effects  are  not  present.  The  third 
term  is,  however,  a  result  of  the  coherent  interference  and  is  respon¬ 
sible  for  forming  the  reconstructed  images.  The  analysis  is  gr'eatly 
simplified  if  it  is  assumed  that  the  developed  photographic  film  is  char¬ 
acterized  by  a  linear  relationship  between  the  exposure  intensity  and 
the  amplitude  transmission  coefficient.  Although  such  a  relationship 
rarely  exists,  the  linear  terms  are  at  least  present  in  the  general  ex- 
•  pansion  relating  the  two.  The  assumption  of  linearity  furthermore 
requires  a  modulation  transfer  function  of  unity,  a  condition  which  is 
never  given  by  photographic  materials.  The  assumption,  therefore, 
represents  a  "best  possible  case"  limit.  The  equations  hereafter  are 
based  on  the  linearity  assumption. 


Expanding  the  sine  function  yields  two  terms,  ip  real  and  ip  virtual, 
the  real  and  virtual  focusing  terms,  respectively,  where 


and 


i  a2  kj 

tf'real  =  | - 

z0  ZJ  (z0  +  Zj) 


— i  a2  kx 

^virtual  =  -  ■ 

z0  zx  (z0  +  Zl) 


exp 


exp 


r  lti- 

L  2zj  (z„  + 

i  kx  z0 
2zx  (z0  +  zx) 


zi) 


■'] 


'•] 


(3) 


(4) 


Only  the  real  image  term  will  be  used  in  the  following  analysis  since  it 
can  be  recorded  without  lenses.  The  background  and  the  virtual 
image  term  can  be  ignored  as  long  as  the  far -field  condition  (Eq.  (1)) 
is  observed. 


The  light  distribution,  0(w)(Fig.  1),  in  a  plane  apart  from  the 
hologram,  formed  by  light  from  the  point  source  s',  which  is  diffracted 
by  the  hologram,  is  given  by  the  Fresnel-Kirchoff  diffraction  integral 
(Ref.  13) 

ipivi)  =  c  /aperture  'Pi  ^  D(r)e  exp[ik2R  (r,  w)]dF  (5) 


where 


c 


(1  —  i  kz  Z2) 
277  Zj2 


cos  9 ' 


-^2 

277Z2 


<6) 
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(where  the  paraxial  approximation  has  been  used)  and 


Z2  =  distance  from  the  hologram  plane  to  the  observation  plane 

d'  =  the  angle  between  the  line  joining  a  point  on  the  hologram 

to  a  point  in  the  observation  plane  and  the  normal  to  the 

hologram  (cos  6'  =  -  =  1) 

R 


^j(r)  =  light  distribution  incident  on  the  hologram 
D(r)  =  amplitude  transmittance  of  the  hologram 


The  light  from  the  point  source  incident  on  the  hologram  is  from  the 
quadratic  approximation 


=  ~  ““T —  exp  (ik2R.2)  exp  ?— ) 

2ffR2  2R2/ 


(7) 


where  R2  is  the  distance  from  the  point,  source  to  the  hologram  plane. 
Expanding  R(r,  w)  in  the  Fresnel -Kirchoff  integral  in  a  binomial  ex¬ 
pansion 


R(r,w) 


+ 


+ 


(8) 


and  substituting  Eqs.  (3),  (7),  and  (8)  into  Eq.  (5),  the  Fresnel -Kirchoff 
equation  yields: 


-2] 

aperture  J 

E  T  Gi  -  T,)'1] 


exp 


k2i 


(2n)2  z  2 R 2 


(9) 


(10) 


(w  and  ?  are  the  cartesian  vectors  in  the  reconstruction  plane  and 
hologram  plane,  respectively. )  Expanding  their  dot  product  in  polar 
coordinates  gives 


w  •  r 


wr  cos  (0  —  9J) 


and  also 


dr  =  rdrd5 
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The  integral  now  becomes,  after  separating  the  integration  over  0 
from  that  over  r, 


.  „  rH  i  /kiar\  -i  ki  z0  2  .  k2  /  1  A  .  . 

Ir  (w,z2)  =  K  i0  Jl  ( -  exp  - - -  r  exp  i  —  — -  +  —  r  dr 

V  zl  /  j_zi(z0  +  zi)J  2  \  R  2  2  2/ 


2 n  I  ik 2  wr  cos  {6  -  d)) 

f  exp  - —I  dd 


z2 


-] 


(11) 


Integration  along  the  radial  vector,  r,  has  an  upper  limit,  H,  corre¬ 
sponding  to  the  limiting  aperture  of  the  hologram. 


The  integral  over  the  angle  6  is  recognized  as  being  the  integral  repre¬ 
sentation  of  the  zeroth-order  Bessel  function,  i.  e. , 


2*7 

J  exp 


[— irw  ,  _  /  k2rw  \ 

2 7T  cos  (  d  =  2nio 


Therefore, 


H  /klar\  /  k2rw  \  f  ->kiz0rZ  “| 

*w’22>  -  &K  ;  Ji  f—j  J.(— )  «* 
ipH.('_L  +  -L')Jdr 

l_2  \r2  z2y  j 


exp  h 


(12) 


(13) 


This  integral,  evaluated,  describes  the  light  distribution  in  a  plane  Z2 
from  the  hologram.  Because  it  is  independent  of  (f> ,  the  bar  has  been 
left  off  the  w,  indicating  that  it  now  is  a  radial  vector. 


Using  Euler's  Formula 

exp  (ig)  =  cos  d  +  i  sin  d 


and  the  fact  that 


I(w)  =  ^(w)  ^*(w) 


(14) 


(15) 


the  intensity  in  the  plane  at  Z2  is 


I(w,Z2) 


(2n)2 


Jl 


(Dr)  J0  (vr)  cos 


+ 


Jl  (Dr)  J0  (vr) 


(16) 
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where 


ki-a 

D  =  — — 

zi 

(17) 

kow 

V  =  — 

(18) 

.'-k,/1  _  >\_ 

\R2  Z2  J  zj(z0  +  zj) 

(19) 

Before  proceeding  further  to  the  case  of  plane  wave  recording  and 
reconstruction,  it  will  be  shown  that  this  equation  indeed  produces  an 
enlarged  image  of  the  particle.  Let  u'  =  0  which  is  the  focusing  condi¬ 
tion  for  the  image  and  assume  the  hologram  contains  all  of  the  diffracted 
light,  i.  e. ,  H  =  <30.  Then  Eq.  (16)  becomes 


I(w,Z2) 


(2ff)2  |K|2  j/“ 3i  (Dr)  J0  (vr)  dr 


(20) 


The  solution  to  this  integral  is  (Ref.  14) 


I(w,z2)  =  (2*)2  |K  2 


for 


for 


for 


lc2w  Ici.a 


< 

z  2 

Z1 

k2W 

kja 

z  2 

Z1 

k  2  w 

,  > 

k  la 

z  2 

Z1 

(21) 


The  edge  of  the  particle  image  is  defined  by 


w 


ki]  Z2 
k2  zi 


(22) 


and  the  particle  is  magnified  by 

ki  Z2 

m  =  T7  TT  (23) 

The  image  location,  Z2,  is  found  from  the  focusing  condition  and 
Eq.  (19)  and  is 

ki  _  kg 
k2zi  k2K  i 
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where  Rj  =  aj  +  zc. 

Thus,  when  the  integral  is  over  infinity  and  the  focusing  condition 
is  applied,  the  magnified  particle  image  reconstructs  at  Z2- 


The  special  case  of  recording  and  reconstructing  the  hologram  with 
plane  waves  can  be  found  by  multiplying  the  spherical  wave  intensity  by 
z04  and  R22  and  taking  the  limit  as  z0  and  R2  go  to  infinity.  The  reason 
for  multiplying  by  zG^  and  R2^  is  seen  by  noting  that  plane  waves  are 

described  by  A0  exp(ikz)  while  spherical  waves  are  described  by  —  ,exP^kr) 

9  A2 

In  the  intensity  equation  these  terms  become  A0^  and  — - ,  respectively, 

^  r>Z 

2  L 

and  the  intensity  equation  must  be  multiplied  by  R2  to  remove  these 
terms  in  the  denominator.  The  term  zQ  is  to  the  fourth  power  since  the 
intensity  is  equal  to  the  square  of  the  integral  containing  the  amplitude 
transmittance  of  the  film,  and  the  transmittance  is  proportional  to  the 
intensity  recorded.  The  intensity  resulting  from  plane  waves  is,  there¬ 
fore, 

Ip.w.  (w,Z2)  =  Lim  R22  zo*  Is.w.  (w,z2) 

R  2  -*  “ 


c2  |j/o  Ji  (Dr)  Jo  (vr)  cos  r2^  dr J 
|^/o  Jl  (Dr)  Jo  (vr)  sin  (~  r2^  drj  | 


(25) 


where 

k22a 

C  2?7z  2 


and 


u 


k2 

z2 


kl 

21 


(26) 


(27) 


Although  this  is  a  particular  condition  of  the  more  general  spherical 
wave  case,  it  has  been  found  that  data  reduction  of  a  volume  of  particles 
is  more  easily  accomplished  with  plane  waves  because  plane  waves  are 
easier  to  reproduce  for  the  reconstruction  and  because  the  magnification 
will  now  always  be  unity  for  any  k2/ki  ratio.  The  focusing  condition 
implies  that 

72  =  +  if  21  (28) 
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The  positive  sign  indicates,  by  convention,  that  the  image  is  located  to 
the  right  of  the  hologram  plane. 


The  integral  in  Eq.  (25)  will  be  known  as  the  focusing  integral.  Its 
solution  for  any  and  all  w  and  Z2  can  be  done  only  by  numerical  integra¬ 
tion.  The  computer  program  which  was  written  to  perform  the  integra¬ 
tion  is  given  in  Appendix  III.  In  order  to  simplify  the  limiting  param¬ 
eters  of  the  hologram,  the  following  special  cases  have  been  examined: 


1.  =  k2  =  k  The  same  wavelength  is  used  in  recording 

and  reconstruction. 


2. 

3. 


Z2  =  zi  +  t]  The  term  rj  defines  planes  to  the  left  and  right 
of  focus  which  occur  when  Z2  =  z.\  (Fig.  1). 


£2Azi 

H  = - —  where  £2  is  the  argument  of  the  first-order  Bessel 

2?ra 


Function  J 


■(“-Sr) 


at  the  upper  limit,  H,  of  the  integral. 


Therefore,  the  limiting  aperture  of  the  hologram  can  be 
defined  in  terms  of  the  number  of  zeros  of  Ji(£2)  recorded. 


4.  The  F-number,  F /#,  of  the  hologram  is  F /#  =  -gg  =  (29) 

This  parameter  is  important  in  that  it  allows  the  hologram 
imaging  properties  to  be  compared  to  a  thin  lens. 


The  integration  method  in  the  program  is  Gaussian.  A  32-  or 
96 -increment  integration  is  available  from  an  input  command,  the 
96 -increment  Gaussian  method  being  more  accurate.  The  number 
of  sine  peaks,  n,  within  the  upper  limit  of  the  integral  should,  for 
the  sake  of  accuracy,  be  less  than  the  number  of  increments,  where 

Input  data  are  the  recording  distance,  wavelength.  A, 


n 


2*2' 


far-field  number,  N,  and  £2,  the  argument  of  Jj(£2),  which  defines  the 
upper  limit  of  integration,  H.  Other  input  data  are  the  minimum  and 
maximum  values  of  rj/a  (rj  normalized  to  the  particle  radius)  and  the 
increments  of  rj/a,  and  the  minimum  and  maximum  values  of  w/a  and 
its  increments.  Output  data  consist  of  t.\,  A,  £2,  N,  a,  F/#,  H,  and 
1(0,  0)  on  the  first  data  page.  The  term  1(0,  0)  is  the  intensity  in  the 
plane  of  focus  at  the  center  of  the  particle.  It  is  the  value  to  which 
the  other  intensities  are  normalized.  The  remaining  data  are 
*norm^  tabulated  for  each  w/a  and  r\j a. 
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2.2  SPECIAL  CASES  OF  THE  FOCUSING  INTEGRAL 

In  order  to  check  the  data  from  the  computer  program  and  also  to 
analyze  the  output  data  and  their  trends,  special  cases  of  the  focusing 
integral  were  evaluated.  These  cases  for  the  most  part  were  chosen 
because  a  closed  form  solution  of  the  focusing  integral  could  be  derived. 

2.2.1  Case  1  -  Ideal  Image  of  a  Perfect  Hologram 


This  case  has  already  been  studied  for  the  spherical  wave  case, 
and  the  results  will  be  repeated  here  for  plane  waves.  In  the  focal 
plane,  at  u  =  0  and  H  =  00 ,  the  particle  image  intensity  is 

for  w  <  a 

for  w  =  a 

for  w  >  a  (30) 

The  ratio  of  the  edge  intensity  to  the  central  intensity  is  seen  to  be  1/4. 
At  the  edge  of  the  particle,  i.  e. ,  the  point  of  discontinuity,  the  ampli¬ 
tude  is  found  by  taking  the  average  of  the  upper  and  lower  values  (1/2). 
The  intensity,  therefore,  is  the  square  of  this  value,  or  1/4  the  inten¬ 
sity  on  the  image  center. 

2.2.2  Case  2-  Intensity  at  the  Image  Center 

All  of  the  computer  data  are  normalized  to  the  intensity  at  the 
center  of  the  particle  image.  For  this,  case,  where  r)  =  w  =  0,  the 
integral  becomes 

I(o,o)  =  -4-  IYji  (Dr)  dr]  -  -4  E  -  Jo  (0)| 

D2  L°  J  D2L  j  (31) 

Thus,  as  the  limiting  aperture  varies,  the  central  intensity  of  the 
particle  image  will  oscillate. 
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2.2.3  Case  3 -Intensity  at  the  Image  Edge 


is  . 


For  a  finite  aperture,  the  intensity  at  the  edge  of  the  particle  image 


I(a,o) 


Jo  (Dr)  Ji  (Dr)  dr 


"il 

r 


(32) 


where  v  =  D  for  w  =  a. 


From  Gradshteyn  and  Ryzhik  (Ref.  14),  the  solution  to  this  integral 


is 


It 


2  r~  l2 


(33) 


Note  that  the  intensity  at  the  edge  of  the  particle  is  a  function  of  the  lim¬ 
iting  aperture  only.  As  the  aperture  becomes  large,  it  approaches  the 
value  given  by  Eq.  (30)  of  Case  1  as  a  limit.  The  intensity  variation 
I(a,  o)  normalized  to  this  limit  is  shown  in  Fig.  2.  The  minima  occur 
when  Ji  (fi)  goes  to  zero,  and  the  peaks  occur  when  J0(£2)  is  zero.  It 
has  been  shown  (Ref.  15)  that 


2  |i  Jr2  (ft)  =  l  -  Jo2  (ft) 


(34) 


Thus,  when  is  a  zero  of  J0(f2),  the  summation  is  a  maximum.  Dif¬ 
ferentiating  with  respect  to  gives 


dO 


[s,jH  - 


2  J0  (ft)  Ji  (ft) 


(35) 


Thus,  the  zeros  of  Ji  (ft)  are  the  minima  of  I.  Figure  3  is  the  intensity, 
I(o,  a),  normalized  to  the  intensity  at  the  center  of  the  particle,  1(0,0). 


2.2.4  Case  4 -Gradient  of  Intensity  at  the  Image  Edge 


The  slope  of  the  intensity  at  the  edge  of  the  in-focus  particle  can 
also  be  derived  as  a  function  of  the  limiting  aperture.  The  intensity  of 
the  particle  in  focus  is  given  by 


”  H  "  ij 

I(w,o)  =  c'2  /  Ji  (Dr)  J0  (vr)  dr 

O 


(36) 
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Taking  the  derivative  of  I(w,o)  with  respect  to  w  gives 

~~  I(o,w)  =  2c2  £/  Jo  (vr)  J 1  (Dr)  drj  ^ — j  £  r  J  i  (Dr)  J0  (vr)  dr 

The  integrals  will  be  evaluated  when  w  =  a,  therefore, 

d  /oc2k  \  H  H 

-  I(o,w)w_.  =  -( -  )  /  Ji  (Dr)  Jo  (Dr)  dr  /  Ji2(Dr)rdr 

dw  m\  zi  I  a  o 


(37) 


(38) 


The  second  integral  has  a  closed  form  solution  (Ref.  15)  given  in 
general  by 


A 

/  z  Jv2  (kz)  dz 
0 


1 

2 


A2  Jv2  (kA)  -  Jv-1  (kX)  Jv+i  (kX) 


(39) 


Therefore, 


H  2 

/  r  Ji2  (Dr)  dr  =  -S-  Jj2  (HD)  -  J0  (HD)  J2  (HD) 

o  2 


Substituting  this  into  Eq.  (38)  and  recognizing  the  first  integral  as 
being  the  same  as  was  solved  in  Case  3,  Eq.  (32),  and  that  HD  =  ft, 
then 

—  I(o,w)w  =  a  =  -  —  1  J2t1  (fi)  (Q)  -  Jo  (D)  J2  («) 

dw  zj  D  n— ° 

=  k12fl2  I  J„  +  1  (D)  J!2  (fl)  -  Jo(Q)  J2(«)  (40) 

n=o 

where 


The  slope  at  the  edge  of  the  particle  is  seen  to  increase  as  the 
aperture  radius  squared.  The  graph  of  this  equation  is  shown  in  Fig.  4. 
The  "staircase"  effect  in  Fig.  4  is  a  result  of  the  Bessel  functions. 

The  increase  in  the  slope  is  zero  for  ^  equal  to  a  zero  of  the  first-order 
Bessel  function.  Normalizing  the  slope  to  the  intensity  1(0,  0)  produces 
the  curve  in  Fig.  5.  Here  it  is  shown  how  the  slope  at  w  =  a  can 
actually  be  less  for  larger  ft  values  than  for  the  smaller  ones.  The 
minima  and  maxima  of  the  curve  occur  at  zeros  of  Jj_.  This  curve  also 
approaches  infinity  as  ft  approaches  <=.  One  other  important  point 
brought  out  in  the  normalization  process  is  that  the  slope  of  this  curve 
is  inversely  proportional  to  the  radius,  a.  Thus,  the  edges  of  the 
smaller  particles  will  be  better  defined  for  a  given  recorded  ft. 
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2.2.5  Case  5- Intensity  along  the  Image  Axis 


As  a  limiting  case  for  the  intensity  away  from  the  plane  of  focus, 
let  H  — ►aD  and  w  =  o,  i.  e. ,  the  intensity  as  a  function  of  rj  along  the 
particle  axis.  This  will  be  called  the  optic  axis  of  the  particle  since 
the  light  distribution  is  symmetrical  about  this  axis  for  each  particle. 
This  is  evident  from  the  fact  that  the  intensity  integral  is  independent 
of  the  polar  angle  $  .  With  these  assumptions,  the  focusing  integral 
becomes 

J  i  (Dr)  drl  + 


I(n,o)  =  c‘ 


i  cos  (t  r2) 


From  Ref.  14  the  closed  form  solutions  of  the  sine  and  cosine 
integrals  are 


{  .MDrHr  =  ±  sin 

/°°cos  /JLr2\  j  (Dr)  dr  =  —  sin2  (— ) 
O  \  2  /  D  \  / 


(42) 


where  >  0  and  D  >  0. 
gives 


I(n,o)  = 


Substituting  these  solutions  into  the  equation 


(43) 


Using  the  trigometric  identity 

2  sin2  —  =  1  -  cos  x 
2 


the  intensity  along  the  optic  axis  of  the  particle  becomes 


I(n,o) 


(44) 


where  D  and  u  have  been  replaced  by  their  definitions  and  p  >  0.  This 
is  a  sinusoidally  varying  function  which  increases  in  frequency  as  r? 
approaches  0.  The  frequency  v  as  a  function  of  rj  is 

77  7.  I 

BN  rj 2  (45) 
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The  last  minima  and  maxima  from  the  focal  plane  occur  when 


l)  mu\ 


z  1 


IN  -  1 


zi 

-IN 


^min 


'\  *1 
8N  -  1  ~  8N 


(46) 


Figure  6  illustrates  the  variation  of  I  along  the  axis  for  N  =  1*  10, 
100.  The  intensity  is  normalized  to  the  intensity  at  1(0,0)  which  is  sim¬ 
ply  a  constant,  for  f2  =  ®.  The  horizontal  axis  is  in  terms  of 

p/a^  to  allow  the  curve  to  be  drawn  for  all  particle  radii  and  recording 
distances.  As  the  particle  image  is  brought  into  focus,  the  intensity 
along  the  axis  oscillates  rapidly  with  z. 

Cases  2,  3,  and  4  apply  equally  well  to  spherical  wave  analysis  since 
the. argument  J^(^)  is  the  same  for  plane  and  spherical  waves.  The  edge 
of  the  image  would  now  correspond  tow  =  ma,  where  the  constant  in  front 
of  the  solutions  to  the  integrals  will  be  the  only  difference.  For  Case  5, 
u  is  replaced  by  u'  for  spherical  waves,  and,  therefore,  the  frequency  of 
oscillation  and  the  points  of  minima  and  maxima  will  be  functions  of  the 
positions  of  the  sources  s  and  s',  and  the  particle  position. 


2.3  DATA  FROM  THE  FOCUSING  INTEGRAL 

The  intensity  distribution  about  the  plane  of  focus  is  a  function  of  a 
number  of  parameters  (e.  g.  ,  particle  size,  recording  distance,  and 
limiting  aperture).  Therefore,  in  order  to  obtain  curves  from  the  com¬ 
puter  data  which  have  meaning  for  the  largest  number  of  cases,  the 
parameters  u  and  v  in  the  intensity  integral  were  chosen  for  the  horizon¬ 
tal  and  vertical  axes,  where 


A  fiz  2  l 

- r  - 

277  a  z\  +  T) 

n  zi 

a  *  z\  +  r\  (47) 


As  is  evident  from  these  equations,  u  and  v  are  not  independent  for 
rj  =  f(zj).  They  are  linearly  related  as 


V 


=  (Jl  —  — - —  u)  m 
2Nfl 


(48) 


where  m  =  w/a.  The  intensity  along  any  line  parallel  to  the  optic  axis 
can,  therefore,  be  constructed  on  a  graph  and  the  intensity  variation 
observed  as  n  varies.  Isophote  diagrams  (i.  e. ,  lines  of  constant  inten¬ 
sity)  are  used  to  display  the  computer  data.  For  a  constant  each 
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isophote  diagram  represents  the  conditions  for  any  particle  radius,  a, 
and  recording  distance,  zj_,  as  well  as  wavelength.  Only  positive  v  is 
shown  in  all  the  diagrams  except  for  one,  since  the  intensity  variation 
for  v  is  symmetric  about  the  optic  axis.  The  diagram  is  also  symmetric 
about  the  v  axis  since  |a  +  ibp  =  |a  -  ibp. 

Some  typical  isophote  diagrams  are  shown  in  Figs.  7  through  11. 
Dashed  lines  on  the  diagrams  indicate  the  intensity  variation  along  a 
line  adjacent  to  the  particle  edge  and  parallel  to  the  optic  axis  as 
illustrated  in  Fig.  7.  This  line  is  found  by  letting  m  =  1  in  Eq.  (48). 

For  small  as  in  Fig.  7,  where  D  =  3.  832,  i.  e. ,  the  first  zero  of 
J0(f2)  /(the  Airy  Disk),  the  intensity  variation  at  w  =  a  is  seen  to  be 
extremely  unsymmetrical  as  predicted  by  Eq.  (4).  More  will  be  said 
about  this  later  when  in-focus  criteria  for  the  image  will  be  defined. 

Figure  8  is  the  isophote  diagram  for  ^  =  5.  0.  The  intensity  varia¬ 
tions  are  shown  to  be  more  complex  about  the  focal  plane.  Figures  9 
and  11  correspond  to  the  second  and  fourth  zeros  of  Ji(^),  respectively. 
Figure  10  illustrates  the  complexity  of  intensity  variations  as  becomes 
larger.  Here  ^  =  11.5.  To  illustrate  the  360-deg  intensity  variations, 
all  four  quadrants  of  the  u  and  v  space  are  shown. 

To  better  illustrate  the  variations  of  intensity  for  various  particle 
radii  and  recording  distance,  graphs  of  I(w,rj)  versus  w/a  and  77/ a  are 
presented.  The  coordinates  w/a  and  rj/a  are  used  to  indicate  the  dis¬ 
tance  in  particle  radii  from  its  center.  Figures  12  to  16  illustrate  the 
intensity  variations  for  typical  particle  radii  and  recording  distances  in 
selected  planes  out  of  focus.  These  would  be  the  intensity  variations 
one  would  observe  over  the  image  of  the  particle  as  planes  of  focus 
were  observed. 

Figures  17  to  ,21  show  the  intensity  variations  along  rj  for  the  optic 
axis  and  the  edge  of  the  particle  as  the  image  is  brought  in  and  out  of 
focus.  These  curves  will  be  used  to  define  the  criterion  for  particle 
size  resolution  and  the  focal  tolerance  (depth-of -field)  of  the  hologram. 


2.4  VALIDITY  OF  SCALAR  DIFFRACTION  THEORY 

The  expressions  derived  so  far  are  all  based  on  a  scalar  diffraction 
theory,  which  is  not  valid  in  general.  Under  certain  conditions  the  vector 
nature  of  the  electromagnetic  wave  must  be  taken  into  account.  This 
occurs,  for  example,  at  small  F /#  ratios.  It  also  occurs  under  some 
polarization  conditions.  Finally  scalar  diffraction  is  not  applicable  in 
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the  proximity  of  a  boundary.  The  complexity  of  the  more  general  vector 
diffraction  theory  encourages  the  use  of  the  scalar  approximation,  which 
should  be  acceptable  for  most  of  the  cases  covered  here. 


SECTION  III 

THEORETICAL  RESOLUTION 


An  accurate  determination  of  size  and  position  of  the  reconstructed 
image  is  important  for  measuring  the  velocity  and  size  distribution  of 
particles  in  a  volume.  From  the  foregoing  data  a  number  of  techniques 
may  be  postulated  by  which  to  determine  the  ultimate  accuracy  of  such 
an  instrument.  The  limits  of  resolution  which  will  be  formulated  will  be 
entirely  theoretical.  It  must  be  remembered,  therefore,  that,  in  prac¬ 
tice,  the  experimental  data  can  be  worse  than  this.  These  will  give  an 
indication  of  the  theoretically  achievable  resolution. 


3.1  PARTICLE  POSITION 


is 


The  smallest  particle  that  can  be  resolved  by  the  in-line  hologram 


amin 


0.6  lAz  i 
H 


(49) 


where  it  was  assumed  that  the  first  zero  of  the  first-order  Bessel  func¬ 
tion  is  recorded.  The  data  from  the  isophote  diagram  for  ^  =  3.  832 
corresponds  to  this  minimum  resolvable  radius. 

The  in-focus  particle  image  will  not,  in  general,  be  the  same  size 
as  the  particle  itself  even  with  plane  wave  recording  and  reconstructing 
beams.  This  is  because  the  edge  of  the  particle  image  has  a  finite  slope 
and,  as  shown  in  Figs.  4  and  5  and  Figs.  12  to  16,  the  intensity  at  the 
edge  of  the  particle  is  not  the  same  for  different  limiting  apertures,  fi. 
Each  particle  will  have  a  different  associated  with  it  according  to  its 
size  and  distance  from  the  hologram.  Thus,  a  criterion  must  be  chosen 
which  will  enable  an  instrument  to  measure  the  particle  size.  The  crite¬ 
rion  which  will  be  used  is  the  intensity  limit  at  the  particle  edge  given 
by  Eq.  (30),  i.  e. ,  Inorm  =  0.  25  where  r?  =  ®.  Figure  3  shows  that  the 
intensity  oscillates  about  this  value  and  approaches  it  as  the  radius  of 
the  limiting  aperture  increases.  The  edge  will  be  assumed  to  lie  at  the 
point  where  the  intensity  is  25  percent  of  that  at  the  center.  The  accuracy 
of  the  measurement,  therefore,  will  depend,  theoretically,  on  how  much 
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of  the  Bessel  function  can  be  recorded,  and  the  maximum  error  will  be 
a  function  of  ft.  Talcing  the  derivative  of  lnorm^a»  °)  with  respect  to  ft 
shows  that  the  extremes  occur  at  the  zeros  of  J ^ (f2 ) .  From  the  isophote 
diagrams.  Figs.  7  to  11,  the  distance  Av,  corresponding  to  the  change 
in  intensity  at  the  known  edge  to  where  it  is  25  percent  of  the  center 
intensity,  was  found.  Using  Eq.  (47)  this  value  can  be  expressed  as  a 
percent  of  particle  radius,  as 

Am  =  A  (v)  "ff  x  100  percent  (50) 

It  is  shown  that  the  error  is  the  same  for  all  particle  radii  at  a  given 
constant  recorded  ft.  Table  I,  Appendix  II,  lists  the  percent  error  for 
typical  ft.  The  errors  designated  as  maximum  correspond  to  the  ft  at 
the  extreme  of  the  curve  in  Fig.  3. 

All  limiting  apertures  above  ft  =  7.  016  have  images  which  contain 
peaks  making  possible  a  second  criterion  for  the  maximum  possible  the¬ 
oretical  error  of  the  images.  These  maxima  occur  just  before  the  edge 
and  exceed  the  intensity  at  the  center  of  the  particle  image.  This  is 
shown  in  Figs.  12  to  16.  These  peaks  can  be  detected  and  the  particle 
size  can  be  found  to  within  this  accuracy.  The  errors  arising  from 
this  measurement  are  given  in  Table  II.  These  values  were  taken  from 
the  isophote  diagrams  for  the  respective  values  of  ft.  Note,  however, 
that  for  ft  =  5  and  3.  832  there  are  no  peaks  involved.  Thus,  the  error 
involved  in  these  cases  is  entirely  a  function  of  the  ability  of  the  detector 
to  measure  the  25 -percent  point.  On  film,  if  the  exposure  level  of  the 
recorded  image  is  small,  only  the  relatively  large  central  peak  of  the 
intensity  will  be  recorded. 

The  particle  position  is  important  since  it  completes  the  data  for  the 
particle  size  distribution  in  the  volume  and  when  the  recording  is  taken 
with  a  double  pulse,  it  gives  the  velocity  of  the  particle.  The  position  of 
the  particle  in  the  plane  of  focus  (i.  e.,  the  p,  q,  plane  of  Fig.  lb)  is 
accurate  to  the  limits  set  by  the  Rayleigh  criterion,  i.  e.  ,  when  the  prin¬ 
cipal  intensity  maxima  of  one  particle  lies  at  the  first  minimum  of  the 
other  the  particles  are  just  resolved.  For  plane  waves,  therefore,  the 
particle  separation,  Sp,  is  (Ref.  4) 

0.6  lAz 

Sp  - - - -  =  2.44  Na  (51) 

However,  the  plane  of  focus  of  the  particle  (occurring  when  n  =  0)  is 
less  accurately  known  than  this.  The  accuracy  limits  set  by  the  hologram 
will,  therefore,  define  the  focal  tolerance,  A 77,  of  the  hologram. 
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To  find  the  At]  resolution  of  the  reconstructed  image,  the  isophote 
diagrams  and  Figs.  17  to  21  will  be  used.  From  these  data  a  number 
of  methods  of  measuring  the  focal  tolerance  (depth  of  field)  of  the 
in-line  hologram  are  evident.  The  first  consists  of  measuring  the  inten¬ 
sity  change  on  the  image  as  it  is  brought  into  and  out  of  focus. 

The  point  that  would  be  observed  could  be  at  the  particle  center, 
its  edge,  or  even  a  slight  distance  away  from  the  edge.  Another  method 
would  be  to  measure  the  slope  at  the  edge  of  the  particle  and  determine 
when  it  becomes  maximum.  Measuring  the  slope  change  would  entail 
an  instrument  continuously  scanning  the  particle  as  it  is  moved  out  of 
focus  and  a  logic  network  to  compare  one  slope  measurement  to  the 
preceding  one.  Conceivably  this  can  be  done  with  a  television  camera 
and  some  logic  circuitry. 


At  the  center  of  the  particle  image  the  criterion  that  a  change  in 
intensity  of  20  percent  can  be  detected  will  be  used.  From  the  definition 
of  the  horizontal  axis,  u,  of  the  isophote  diagrams,  the  change  in  r)  is 


At)  =  ± 


2nl 

A 


(52) 


The  value  of 


—jean  be  taken  from  the  isophote  diagrams. 


The  focal 


tolerance,  Arj,  therefore,  will  apply  for  all  particle  diffraction  patterns 
having  the  same  upper  limit  of  the  Bessel  function.  Equation  (52)  can 
be  rewritten  as 


AtI  =  ±  kholo  (ft)  77  (53) 

where 

kholo  (ft)  “  —7-  A  ( — ) 
n2  \n/ 


The  values  of  Kholo(^)  f°r  each  S7  are  given  in  Table  III  along  with 

the  corresponding  from  the  isophote  diagrams.  Notice  that  the 

accuracy  improves  as  ft  increases.  For  ft  =  7.  016  and  X  =  6328  A,  the 
focal  tolerance  of  a  hologram  of  a  50 -Mm  particle  is  ±635  Mm.  However, 
for  a  10-Atm  particle  the  focal  tolerance  reduces  to  ±25.4  Mm. 


All 


— j,  except  when  ft 


3.832,  were  measured  at  a  20-percent  increase 


in  the  intensity. 
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An  ideal  thin  lens  imaging  a  point  source  has  a  focal  tolerance 
given  by  (Ref.  13) 


Ajj  =  ±  3.2  —  (F/ #) 2 

77 


(54) 


where  F /#  is  the  F  number  of  the  lens.  To  compare  this  to  the  holo¬ 
gram,  the  F/#  of  the  thin  lens  will  be  assumed  to  be  the  same  as  that 
of  the  hologram.  Substituting  the  definition  of  the  F/§  for  the  hologram 
(Eq.  (29))  into  this  equation  gives 


A  T)  =  ± 


20.1 

fl2 


-r  ■  kT-L-<a)  -T 


(55) 


The  thin-lens  focal  tolerance  is  of  the  same  form  as  that  obtained  from 
the  hologram  since  the  same  coordinate  relationship  to  17  was  used  in 
both  cases.  The  value  of  Kji  l.  (^)  f°r  the  thin  lens  is  given  in  Table 
III.  It  is  shown  that  the  criterion  chosen  for  the  hologram  yields  results 
for  the  A 17  position  of  the  particle  which  are  comparable  to  a  thin  lens. 

The  intensity  change  at  the  center  of  the  image  was  shown  in 
Section  2.  2.  5,  where  =  “  is  an  oscillating  function  of  17.  For  this 
limiting  case,  a  detector  would  not  need  to  recognize  an  intensity  change 
of  a  given  amount.  Rather,  the  spatial  frequency  of  the  intensity  could 
be  observed  until  the  cutoff  frequency  of  the  detector  or  its  associated 
electronics  was  reached.  This  method  could  conceivably  further  in¬ 
crease  the  resolution.  For  the  small  limiting  apertures,  the  errors 
associated  with  measuring  the  position  of  this  extrema  at  the  center  of 
the  particle  are  large.  The  values  of  Kh0i0(£2)  for  this  criterion  are 
also  given  in  Table  III. 

Observing  the  intensity  change  at  the  edge  of  the  particle  has  its 
inherent  difficulties,  the  first  of  which  is  finding  the  true  particle  edge. 
As  was  shown  previously,  this  value  can  have  a  theoretical  error  of  as 
much  as  20.  5  percent.  Also,  as  indicated  by  the  dashed  lines  signifying 
a  line  tangent  to  the  edge  of  the  particle  and  parallel  to  the  optic  axis, 
the  intensity  is  not  symmetric  on  either  side  of  focus,  except  for  large  £2 
and  N.  A  worse-case  condition  was  shown  in  Fig.  12  where  £2=3.  832 
and  N  =  5.  The  intensity  past  the  point  of  focus  is  much  greater  than  at 
the  focal  point.  This  phenomenon  is  observed  upon  reconstruction  as  an 
intensity  increase  in  the  particle  image  past  the  focal  plane.  If  a  min¬ 
imum  detectable  intensity  change  of  20  percent  of  that  in  the  plane  of 
focus  is  used,  then  Fig.  11,  where  £2  =  13.  324,  indicates  that  this  value 
will  not  be  reached  until  u  >  9tt  which  is  at  least  nine  times  poorer  than 
the  resolution  associated  with  measuring  the  center  intensity. 
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The  third  point  under  consideration  as  an  indication  of  the  Az 
resolution  is  at  w  =  a+  where  a+  is  some  distance  beyond  the  edge  of  the 
image.  When  the  image  is  observed  with  the  eye,  out-of -focus  planes 
are  seen  as  containing  an  Enlarged  image  with  fuzzy  edges.  This 
corresponds  to  the  intensity  beyond  the  particle  edge  increasing.  Here, 
as  in  the  case  of  observing  the  edge  of  the  particle,  difficulties  arise  in 
finding  the  edge  without  a  prior  knowledge  of  the  particle  size. 

Criteria  for  the  particle  size  and  position  have  been  defined  in  terms 
of  an  intensity  change.  The  theoretical  resolution  of  the  particle  size  is 
seen  to  be  better  than  20.  5  percent  for  a  worse-case  condition,  i.  e.  , 
only  the  first  zero  of  the  first  Bessel  function  is  recorded.  Equation  (50) 
shows  that  this  resolution  is  only  a  function  of  the  limiting  aperture, 

On  the  other  hand,  the  depth  of  field  of  the  hologram  is  approximately 
the  same  as  a  thin  lens.  Equation  (52)  shows  that  the  focal  tolerance, 

Arj,  is  a  function  of  the  limiting  aperture  and  the  particle  radius.  This 
can  be  related  back  to  the  F/#  of  the  hologram  to  give 

Arj  =  +  2A  <F/.)2  A  (-£-)  (56) 

The  smaller  the  F/§  is  the  better  will  be  the  focal  tolerance  and,  there¬ 
fore,  Zj  and  the  separation  between  two  image  planes  can  be  more 
accurately  found. 


3.2  LIMITING  PARAMETERS  OF  THE  FILM 

/ 

Up  to  this  point  nothing  has  been  said  about  the  causes  of  the  limiting 
aperture  other  than  the  fact  that  there  is  one  and  that  it  keeps  the  recon¬ 
structed  image  plane  and  the  image  from  being  exactly  known.  The  film 
can  limit  the  aperture  of  the  recorded  diffraction  pattern  in  two  ways. 
First,  because  of  the  linearly  increasing  frequency  of  the  sine  term  in 
the  recorded  intensity  at  some  point  on  the  film,  (Eq.  (2))  the  frequency 
will  exceed  that  which  can  be  recorded.  Secondly,  inherent  on  the  film 
is  its  grain  structure  which  introduces  noise  into  the  recorded  signal. 

The  decreasing  peak  amplitude  of  the  Bessel  function  in  the  interference 
term  in  Eq.  (2)  and  the  1/r  dependence  of  the  intensity  causes  a  reduction 
in  signal  amplitude  as  r  increases.  At  some  point  on  the  film,  rmax, 
the  signal  will  be  engulfed  in  noise  and  will  be  ineffective. 

The  cutoff  frequency  of  the  film  is  related  to  the  recording  param¬ 
eters  through  the  sine  term.  The  derivative  of  its  argument  with  respect 
to  r  is  the  frequency,  v ,  at  a  given  r,  i.  e.  (assuming  an  M.  T.  F.  of 
unity), 

v  =  ~tr  (57) 
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The  radius,  rmax>  corresponding  to  the  maximum  cutoff  frequency  of 
the  film  is,  therefore. 


rmax 


kAzi  =  4Na2v 


(58) 


where  z  = 


N4a" 


was  used.  This  is  the  radius  of  the  limiting  aperture 


set  by  the  flLm  cutoff  frequency.  The  maximum  argument,  ^max> 
the  first-order  Bessel  function  is 


of 


fin 


=  2na  vn 


(59) 


As  an  example,  assume  Kodak®  SO- 243  film  is  used  as  the  recording 
medium;  it  has  a  cutoff  frequency  of  about  200  lines/mm.  For  a  50-Mm 
particle,  therefore,  the  film  will  record  the  Bessel  function  until  its 
argument  is  20?r.  The  size  of  this  aperture  is  2N  mm,  N  being  the  num¬ 
ber  of  far  fields  the  film  is  from  the  particle.  Table  IV  shows  rmax 
and  9max  for  a  range  of  particle  sizes.  The  values  of  rmax  below  the 
heavy  line  will  not  be  limited  by  the  cutoff  frequency  of  a  4-  by  5-in. 
sheet  of  film,  but  rather  its  size. 


The  minimum  F/#  the  film  may  have  is  found  from  Eq.  (29)  and 
Eq.  (59)  and  is  given  by 


F/»n 


(60) 


The  F /#  is  seen  to  be  as  good  as  the  film  cutoff  frequency.  For  SO- 243 
film  and  X  =  6328  A  the  smallest  F/#  attainable  by  this  film  is  3.  95. 

The  Az  resolution  was  shown  previously  to  be  (Eq.  (52)) 


Substituting  Eq.  (59)  for  gives 

A^/min  =  ± 


2it\ 


(61) 


The  resolution  of  the  particle  position  is  seen  to  depend  on  the 
maximum  recordable  frequency.  For  SO-243,  At?  =  ±6.3  x  10"®  4% 


Thus,  this  relatively  low  resolution  film  is  capable  of  resolving  the  focal 
plane  to  better  than  ±10  Pm,  assuming,  that  is,  that  the  cutoff  frequency 
is  the  major  limiting  parameter. 
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From  Eq.  (49)  it  is  seen  that  the  minimum  resolvable  particle  radius 
can  also  be  expressed  in  terms  of  the  cutoff  frequency  as 

0.61 

ami„  «  - -  (62) 

•'m  ik 

From  this  equation,  SO-243  is  capable  of  resolving  particle  radii  as 
small  as  3.  15  Mm. 

The  resolution  capabilities  of  the  hologram  for  small  particles 
appear  from  the  foregoing  criteria  specifications  to  be  extremely  good 
for  a  relatively  low  frequency  film.  However,  as  stated  before,  the 
cutoff  frequency  is  not  the  only  limiting  parameter  of  the  film.  The 
film  grain  noise  which  is  present  on  the  recorded  signal  is  also  a 
limiting  factor. 

The  film  noise  is  measured  in  rms  granularity  and  is  defined  as 
(Ref.  16) 


(63) 


where  ADj  is  the  deviation  of  n  density  readings.  In  practice  the  rms 
granularity  is  measured  by  scanning  an  aperture  of  area  A  over  the  film. 
The  granularity  is  given  to  a  good  approximation  as  (Ref.  16) 


.(D)  . 

v'a 


(64) 


where  D  is  the  average  density  on  the  film  and  G,  which  is  the  Selwyn 
granularity,  is  G  =  2  yj 0 .  43tt  d,  where  d  is  the  mean  diameter  of  a  devel¬ 
oped  grain.  The  granularity  is  seen  to  have  a  reciprocal  relationship 
to  the  scanning  aperture. 


The  signal -to-noise  ratio  (SNR)  of  the  film  is 


SNR 


Db  -  Ds 

<7(D') 


(65) 


where  Dg  is  the  average  density  of  the  background  and  Dg  is  the  average 
density  of  the  signal.  The  normalized  intensity  incident  on  the  film  is 
given  by 


I 


norm 


1 


n 

N 


Ji  <Q) 

n 


(66) 
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where  plane  wave  recording  is  assumed  and  the  third  term  has  been 


dropped  since  it  is  negligible.  Since  ^ 


n  Jl<n> 


modulates  the  sine  term. 


n  n 

the  peak  magnitude  of  this  term  will  only  be  considered.  The  signal 
density  recorded  by  the  film  is  (Ref.  17) 


D  =  y  Log  10  (It)  +  Dc 


(67) 


where  7  is  the  slope  of  the  characteristic  curve  for  the  film,  t  is  the 
exposure  time  and  D0  is  the  density  where  the  linear  portion  of  the 
curve  intersects  the  vertical  axis.  Grouping  the  constants  together 
gives 


D  =  y  Logio  I  +  K 


(68) 


where 


K  =  y  Logio  t  +  D„ 

The  intensity,  I,  can  be  decomposed  into  a  signal  Ig  and  background, 
Ig,  and  D  becomes 


D  =  y  Logio 


+  K' 


(69) 


Here  K'  =  K  +  7  Logio  I3  and  the  signal  intensity  has  been  normalized 
to  the  background.  Since  Ib  >  IS  the  logarithm  can  be  expanded  as 


Logio  (1  +  *) 


1 

2.303 


1 

2.303 


'  X 


Therefore,  the  density  on  the  film  resulting  from  the  plane  waves  and 
diffracted  light  off  the  particle  is 


D  =  - 


2.303N 


Jl  (A) 

n 


+  K' 


(70) 


K'  represents  the  average  background  density  Dg.  Therefore,  Eq.  (65) 
becomes 


SNR 


where  Eq.  (64)  was  used. 


yn  j  1  (fi)  1 

2.303N  fi  (7(D) 

yn  \/T  'Ji(fl) 

2.303NG  VD" 


(71) 

(72) 
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It  is  shown  that  as  the  area  of  the  scanning  aperture  increases  or 
equivalently  when  the  area  of  the  recorded  signal  density  is  large  the 
SNR  is  maximized. 


A  half  period  of  the  sine  function  can  be  assumed  to  be  recorded  on 
the  film  as  an  annular  region  since  the  diffraction  is  symmetrical  about 
the  optic  axis  of  the  particle,  i.  e. ,  its  center.  This  is  illustrated  in 
Pig.  22.  The  area  of  this  region  is 

A  =  77  (r22  -  ri2)  (73) 


where  r^  and  r2  are  the  radii  of  the  inner  and  outer'regions,  respectively. 
The  terms  r\  and  r£  correspond  to  consecutive  zeros  of  the  sine  term. 
Since  T2  =  ri  +  Ar,  Ar  can  be  found  from  this  term  and  is 


ffzii 

kri 


(74) 


where  it  is  assumed  that  Ar  <  ri.  The  area  of  the  annulus  is  now 


A  =  ffAzj 


(75) 


This  is  seen  to  be  independent  of  r  for  a  constant  just  as  the  half- 

N4a^ 

period  zones  of  a  Fresnel  lens  (Ref.  18).  Using  this  area  and  ^ — , 

Eq.  (72)  becomes 


2y7r\fn  a  J  i  (ft) 

SNR  =  -  -  - 

2.303  gVD7Ic  Vn  ft 


(76) 


In  terms  of  the  rms  granularity  of  the  film,  the  average  density  at 
which  it  was  measured,  /Bscan,  and  the  diameter  of  the  scanning  aper¬ 
ture,  daper,  the  SNR  becomes 

4  yn  VDscan  a  Jl(ft) 

SNR  =  -  ...  ■  — —  — “ —  (nn\ 

2.303(7  (D)<l„p  er  v'Drec  VN  ft 

This  equation  and  Eqs.  (71)  and  (72)  represent  the  SNR  of  the  film 
for  the  diffraction  pattern  recorded  on  an  average  background  density 
\  Drec  and  with  a  film  gamma  of  y.  The  quantity  Ji(f2)/f2  is  shown  as  a 
function  of  £2  in  Fig.  23.  The  peak  amplitudes  of  the  first-order  Bessel 
function  were  used  for  this  figure.  It  is  assumed  that  the  signal  is  com¬ 
pletely  lost  in  the  noise  when  the  peak  of  Ji($7)  is  below  the  noise  level. 


For  a  given  SNR  and  granularity  the  limiting  aperture  radius  can 
be  found.  Table  V  gives  the  maximum  Bessel  function  argument,  ^max* 
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recorded  and  the  limiting  aperture  size  of  the  diffraction  patterns  for 
typical  particles  and  far-field  numbers  with  SNR  =  1.  These  are 
assumed  to  be  recorded  on  SO-243  film  which  has  an  rms  granularity 
of  7.  4  x  10“3  measured  with  a  48-Mm-diam  aperture  and  a  mean  den¬ 
sity  of  1  (Ref.  19).  It  has  been  assumed  that  Drec  =  Dscan  and  y  =  2.2. 

The  larger  particles  are  shown  to  be  limited  not  by  the  film  grain  noise 
but  by  either  the  film  size  or  the  size  of  the  beam  illuminating  the 
particle  field.  Maximum  radii  below  the  dashed  line  in  the  table  are 
larger  than  a  4-  by  5 -in.  sheet  of  film. 

Comparing  Tables  IV  and  V,  it  is  seen  that  the  size  of  the  hologram 
on  the  film  can  be  limited  by  either  the  film  cutoff  frequency  or  noise 
depending  upon  the  far-field  number  and  the  particle  size.  In  Table  V 
all  values  of  9max  below  the  solid  line  are  smaller  than  the  maximum 
£2  of  the  cutoff  frequency  for  that  particle  radius.  Thus,  above  this  line 
the  hologram  size  is  limited  by  the  cutoff  frequency. 

Besides  the  effects  of  the  film  grain  and  cutoff  frequency,  system 
noise  is  also  responsible  for  limiting  the  hologram  size.  This  noise 
can  include  dust  and  dirt  in  the  system,  an  inhomogeneous  wavefront 
illuminating  the  particle  volume,  and  improper  film  developing  techniques. 
These  for  the  most  part,  however,  are  not  predictable. 


SECTION  IV 

SUMMARY  AND  CONCLUSIONS 


From  the  intensity  distribution  recorded  in  an  in-line  hologram, 
the  integral  describing  the  intensity  distribution  about  the  reconstructed 
real  image  was  found.  Because  of  its  complexity,  only  after  special 
conditions  were  imposed  upon  this  integral  could  it  be  solved  in  closed 
form.  In  general,  numerical  integration  was  necessary  for  its  solution. 
The  data  from  the  integration  was  put  in  the  form  of  isophote  diagrams 
with  hybrid  u-v  coordinates  corresponding  to  constants  within  the  argu¬ 
ments  of  the  functions  integrated.  These  diagrams  which  illustrate  the 
lines  of  constant  intensity  about  the  particle  image  apply  to  all  cases 
of  particle  radii  and  recording  distances.  The  limiting  aperture  of  the 
hologram,  ft,  or  more  correctly,  the  upper  limit  of  the  argument  of  the 
first-order  Bessel  function,  was  the  parameter  which  was  different  in 
each  of  the  diagrams. 

Resolution  of  the  in-line  hologram  was  found  from  the  isophote 
diagrams  by  specifying  certain  criteria  which  define  the  particle  image 
edge  and  its  in-focus  limits.  Because  of  the  finite  slope  of  the  recon¬ 
structed  image  edge,  two  criteria  were  presented  for  consideration. 
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The  first  took  advantage  of  the  fact  that  the  squared  average  of  the  light 
amplitude  at  the  edge  of  an  image  formed  with  all  the  light  diffracted 
from  the  particle  was  25  percent  of  the  particle  intensity.  The  inherent 
error  in  choosing  this  point  was  seen  to  be  maximum  when  the  aperture 
cutoff  was  at  a  zero  of  the  first-order  Bessel  function.  (The  error  was 
zero  when  the  cutoff  was  a  zero  of  the  zeroth-order  Bessel  function. ) 
Recording  only  the  essential  Airy  Disk  of  Ji(ft),  i.  e.  ,  its  first  zero, 
yields  an  error  of  20.  5  percent,  the  largest  possible  theoretically.  As 
increases  this  maximum  error  reduces  to  zero.  The  second  crite¬ 
rion  was  formulated  on  the  observation  that  the  particle  image  contained 
peaks  in  intensity  preceding  the  edge.  If  the  last  peak  before  the  edge 
could  only  be  discerned  from  the  background,  the  error  in  the  particle 
radius  measurement  would  then  be  45  percent  for  ^  =  7.  016  and  20  per¬ 
cent  for  =  13.324.  As  increases  the  accuracy  also  gets  better. 

This  criterion  could  not  be  used  for  =  3.  832  and  =  5  since  their 
images  contained  only  a  central  peak. 

The  planar  resolution  containing  the  in-focus  image  ideally  should 
be  infinitesimally  small.  This  would  allow  the  image  to  be  located  at  a 
precise  distance,  z^,  from  the  hologram  and  the  distance  between  two 
images,  Az,  could  accurately  be  found  for  velocity  data.  However, 
because  of  the  limited  hologram  size,  the  image  contains  in-focus  char¬ 
acteristics  over  a  finite  distance,  A p.  Minimum  changes  in  these  char¬ 
acteristics,  therefore,  can  be  used  to  define  the  focal  tolerance  or 
depth -of -field  of  the  image.  As  the  image  is  moved  through  focus,  the 
data  showed  that  the  intensity  and  the  slope  at  the  edge  of  the  image 
change.  The  center  of  the  image,  the  particle  edge,  and  a  point  away 
from  the  edge  could  be  observed  as  to  their  intensity  variations. 

Because  of  the  error  inherent  in  locating  the  edge  and  determining  a 
fixed  point  away  from  it  for  all  particles  and  recording  distances,  the 
intensity  at  these  points  and  the  slope  change  were  not  considered. 

At  the  center  of  the  image  the  criterion  that  a  minimum  intensity 
change  of  20  percent  could  be  detected  was  chosen.  The  data  from  the 
isophote  diagrams  showed  that  with  this  criterion  the  hologram  has  the 
same  focal  tolerance,  Arj,  as  a  thin  lens,  and  A n  is  a  function  of  the 
F/#  of  the  hologram.  In  terms  of  the  particle  radius 

Ar?  =  ±KHoio<n> 

where  K(fi)  is  a  constant  for  each  maximum  Ji(fl)  argument,  and  its 
value  is  given  in  Table  III.  For  example,  the  focal  plane  of  a  20 -Pm 
particle  image  can  be  found  to  within  100  pm  if  only  the  second  zero  of 
is  recorded. 
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The  intensity  at  the  center  of  the  image  was  found  to  oscillate  as  it 
is  moved  out  of  focus.  If  the  last  peak  before  the  focal  plane  is  used  as 
the  in-focus  criterion,  the  focal  tolerance  is  approximately  double  that 

found  by  the  20  percent  intensity  change  criterion.  K  . 

J  J  &  hologram  first 

,  in  Table  III  is  seen  in  most  cases  to  be  twice  K,  (fi)... T. 

extrema  holo  20%I 

The  limited  effective  size  of  the  hologram  was  shown  to  be  the 
factor  reducing  the  accuracy  of  the  measurements  of  the  reconstructed 
image  size  and  z  location.  The  parameters  which  cause  this  limiting 
aperture  effect  on  the  film  are  the  cutoff  frequency  of  the  film,  its  grain 
noise,  system  noise  (such  as  unwanted  diffraction  patterns  from  dust  and 
other  particles),  film  size,  and  reference  beam  size.  Diffraction  pat¬ 
terns  of  small  particles  (<50Mm)  recorded  in  small  far-field  distances 
were  found  to  be  limited  for  Kodak  SO-243  by  the  cutoff  frequency.  The 
medium  range  particles  (50  Mm  <  a  <  100 /im)  diffraction  patterns  are 
limited  by  the  loss  of  signal  in  the  grain  noise,  and  the  larger  particle 
(a  >  100 Mm)  patterns  are  limited  by  either  the  film  size  or  the  reference 
beam  size.  System  noise  will  reduce  these  limits  even  more. 

The  focal  tolerance,  Arj,  and  the  minimum  resolvable  particle  radius, 
ajnin,  were  expressed  in  terms  of  the  cutoff  frequency.  This  parameter 
could  be  used  as  a  rule-of-thumb  indicator  for  film  quality  since  both  of 
the  limiting  parameters  of  the  film  are  related  to  a  common  source,  the 
film  grain.  Kodak  SO-243  film  with  a  conservative  cutoff  frequency  of 
200  lines/mm  is  capable  of  resolving  a  3 -Mm  particle  to  within  ±10  Mm. 

Future  work  will  include  experimental  verification  of  the  theoretical 
results  obtained  here.  Low,  medium,  and  high  granularity  film  will  be 
used.  The  results  obtained  for  the  SNR  will  be  compared  to  the  recon¬ 
structed  image  and  the  measured  film  noise.  The  intensity  changes 
about  the  particle  image  will  be  verified.  Criteria  for  clean  recon¬ 
structed  images  will  be  established  to  check  the  correctness  of  the 
theoretical  criteria  chosen. 

Theoretically,  future  work  can  continue  in  improving  the  measure¬ 
ment  of  image  separation  by  studying  interference  effects  between  two 
images  separated  by  a  distance  less  than  their  diameters.  The  theory 
can  also  be  extended  to  spherical  waves  and  the  hologram  resolution  as 
a  function  of  the  total  magnification.  This  could  be  compared  to  that 
obtainable  from  thin-lens  magnification  so  that  an  engineering  criterion 
can  be  evolved  to  pick  between  them. 
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Fig.  5  Normalized  Slope  at  the  Particle  Image  Edge 
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Fig.  12  Intensify  Variation  Perpendicular  to  the  Optic  Axis  for  Various  Values  of  ir 
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Fig.  13  Intensity  Variation  Perpendicular  to  the  Optic  Axis  for  various  values  of  n 
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Fig.  14  Intensity  Variation  Perpendicular  to  the  Optic  Axis  for  various  values  of  n 
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Fig.  16  Intensity  Variation  Perpendicular  to  the  Optic  Axis  for  Various  Values  of  n 
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Fig.  17  Intensity  Variation  at  the  Center  and  Edge  of  the  Particle  Image  as  it  is  Brought  into  and  out  of  Focus 
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Fig.  18  Intensity  Variation  at  the  Center  and  Edge  of  the  Particle  Image  as  it  is  Brought  into  and  out  of  Focus 
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Fig.  19  Intensity  Variation  at  the  Center  and  Edge  of  the  Particle  Image  as  it  is  Brought  into  and  aut  of  Facus 
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Fig.  20  Intensity  Variation  at  the  Center  and  Edge  of  the  Particle  Image  as  it  is  Brought  into  and  out  of  Focus 
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Fig.  21  Intensity  Variation  at  the  Center  and  Edge  of  the  Particle  Image 
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Fig.  22  The  Diffraction  Pattern  Resulting  from  Plane  Wave  Illumination  of  a  Circular, 
Opaque  Particle 
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TABLE  I 

THEORETICAL  ERROR  IK  THE  PARTICLE  RADIUS,  WHEN  THE 
INTENSITY  IS  25 %  OF  THAT  AT  THE  CENTER  OF  THE  IMAGE 


Q 

%  Error 

3.832 

±22.0  %  Max 

5.0 

±8.0  % 

7.016 

±6.0  %  Max 

11.5 

±1.00% 

13.324 

±3.5  %  Max 

TABLE  II 

MAXIMUM  THEORETICAL  ERROR  IN  THE  PARTICLE  RADIUS 
MEASUREMENT  ARISING  FROM  MEASURING  THE 
OUTERMOST  PEAK  ON  THE  IMAGE 


u 

Q 

Maximum  %  Error 

7.016 

47.0% 

11.5 

28. 0% 

13. 324 

19.8% 
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TABLE  III 

FOCAL  TOLERANCE  CONSTANT,  FOR  A  HOLOGRAM  AND  A  THIN-LENS 
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TABLE  IV 

MAXIMUM  LIMITING  APERTURE  RADII  AND  MAXIMUM  ARGUMENTS 
OF  J^D)  FOR  umax  =  200  LINES/MM 


rmax 


a(|im) 

fynax 

\a 

N\ 

10 

50Mm 

100  Mm 

250  Mm 

10 

12.6 

1 

80  Mm 

2  mm 

8  mm 

5  cm 

50 

62.8 

5 

400  Mm 

1  cm 

4  cm 

25  cm 

100 

125.7 

10 

800  Mm 

2  cm 

8  cm 

50  cm 

250 

314.2 

20 

1.6  mm 

4  cm 

16  cm 

1  m 

50 

4  mm 

10  cm 

40  cm 

— 

80 

6.4  mm 

16  cm 

64  cm 

— 

100 

8  mm 

20  cm 

80  cm 

— 
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FORTRAN  It/  G  LEVEL  1.  NQO  3  MAIN  DATE  =  70019  21/06/35  PAGE  0001 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

'0018 

0019 

0020 
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0026 

0027 

0029 

0029 
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0031 
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0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
'0041' 
0042 
0043 
0044 
~  0045 
0046 


IMPLICIT  REAL *8  IA-Ht0-ZI 

REAL*4  XMIN.XMAX.OELTAX, YMIN, YM AX, DELTA V 

COMN0N/KOUNT/  NE.NW 

COMMON/ I NPUT/  XHI N , XMAX , OELTAX, YMIN, YNAX, OELTAY, INDMX, INO 
COMMON  /ALL/  W,O.BL,FK,Z 
COMMON  /ONETA/  ETA 
I C0UNT=0 

395  CALL  ERRSETI207,Z56,5,1I 

RFAO  15,10001  Z, ALAM.GAN.CAPN 
REA0I5, 1000  INN,  MM  AX,  NOEL,  RLC1WER 
REAOI5,1000IETAM,ETAX,E TAD, CODE 
PI-  3.14159765359 
A*  ALAN  4  Z 
A*  A /  (4.0  *  CAPNI 
A*  OSQRTIAI 

WRITEI6, 10201  Z, ALAM.GAN.CAPN 
1020  FORMAT  I 1H1 , ■  Z  =  '.EZO.e,'  LAMBDA-  ',E20.8,/, 
i*  gamma-  '.Ezo.a,*  n=  »,F?n.B)  - 
BL-1.0 


FK-  2.*Pl/ALAM 

FETA  =  ETAM  4  A 

FINFTA-  FTAX*  A 

OELETA  =  ETAD  •  A 

FW  =  NM  *  A 

FINN  =  UMAX  4  A 

OELW  =  NOEL  *  A 

H-IGAH/PI I40S0RTI Z*ALAM*CAPNI 

FN3=  PI/12. 0*GAMI  4  DSORTI  Z / I CAPN*ALAM) I 

WRITE  16, 10211 A,FNO, RLOWER.H 

1021  FORMATUHO,*  A=  *.£20.8,'  F*«  ',620. 8,'  LOWER-  ’,E20.8,*  H=  •, 
4F20.81 

TERM  =  ALAM4Z4  CAPN 
C  ■  2.0  40SQRTITERMI 

0  =  PI  /  OSQRTI TERMI 
FIE=  I ETAX-FTAM) /  ETAO 
F IW-  (WMAX-WM1/W0EL 
NW  =  FIN  41.0 
NE  =  FIE  *1.0 
XMIN-  ETAM 
XMAX-  ETAX 
"OELTAX-  ETAO 
YMIN-  WM 
YMAX-  WMAX 
OELTAY  WO  EL 

eta-o.6  . 

W-0.0 

IFIC0DE.EQ.96. IGO  TO  69 
_  CALL  GSS321 1  ,RLOWF.R,H,E  1 ) 


|  / / AKDO 1549  JOB  (ARO, 

* 

// 

SRSO0452 ,01, BC 59161,'  B322  P/0  020' , M SGLE VEL- 1 , 

* 

// 

CLASS-C 

// 

FXFC 

FIGLNKG0,MCHL1ST=LIST 

XX 

PROC 

BLKF  AC-1'596,  0ATABLK  =  400,PLIITLI  B-PL0T835 , MCHL  I  ST-NOL  1ST 

XXFTG 

FXEC 

PGM-IEYFORT.PARM-’ SOURCF .MAP , EMCHL I  ST • 

XX5YSLIN  DO  UNIT-WORK, 0 I SP  =  ( , P AS S 1 , SPACE - 1 C YL, 13, 11  1, 

XX 

DCB=IRECFM=FB,LRECL=B0,DLKSIZE-32O0) 

XXSY5PRI NT  DO 

S YSOUT-A, UNI T -PRINT, SPACE-I 1200, I 700 , 100 ), RLSE ,, ROUND ) , 

* 

XX 

nCB-lRECFM=FBA,LRECL=I20,BLKSIZE-I200) 

//SYSIN 

OD  4 

IEF236I 

ALLOC. 

FOR  ARD01549  FTG 

IEF237I 

5YSLIN 

ON  130 

IEF237I 

5Y5PRINT  ON  136 

I EF2  37  J 

SYSIN 

ON  131 
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CALL  GSS32I 2.RL0WER . H,E 2  I 
GO  TO  6464 

64  CALL  GSS46I 1, KLOWFR , H,F 1 1 
CALL  GSS96(2.RL3WER,H,E2I 
6464  CONTI NUF 

ETAN»  C**2  *  I E 1**2  +  E2**2I 
WRITF  (6,10101  ETAN 
1010  FORMAT  ( IHO  t  *  1(0,01  =  '.E20.8I 
W  =  rw  -  OELW 
ETA=FETA-  OELETA 
DO  1  1-1, NU 

W*  M  +  OELW 
OD  1  J«1,NE 
ETA  *  FT  A  +DELETA 

IF  ( J.EQ. l.OR.J.rO. 26. OR. J. E0.51.0R. J. EO. 761  CALL  FLIP 

IF(CnOF.e0.56.IG()  TO  46 

CALL  GSS32I  l .RLOWFR , H, E 1 1 

CALL  CSS32(2.RLOWER.H,E2l 

GO  TO  9646 

96  CALL  GSS96( 1, RLOWFR, H,E1) 

CALL  G$S96(2,RL0HFR,H,E2I 
9696  CONTINUE 

ETA  I  =  C**2  *  (El **2  +  E2**2 I 

VNORN—  ETAI/ETAN 
THEW=  WM»  ( ( I- 1 1 *HDE  LI 
THEET  A=  ETAM  *■  C I J- 1 1 *E  T AD  I 
WRITE (6, 10021 THEW, THEETA , VNDRM 
PHI  =  [)ATAN2(E2,E1  I 
AMAG=0S3RT(E1**7»F2**2I 
WRITE (6, 20001 AMAG, PH I 

23  00  FORMAT! 25X, *MAGN ITUDE= ’ ,E 1 5 . B, 5X, • ANGLE ( RADI  ANSI  =  ' ,E 15 . 81 
IF  (J.EO.NEI  ETA*  FET A-DFLI  T A 
l  CONTINUE 

1000  FORMAT  (4E15.01 
1002  FORMAT  {20X.3E20.8I 
I  COUNT* I COUNT* 1 
I F ( IC0UNT.EQ.1IWRI TF(6, II  111  I 
lllll  FORMAT!//,'  RUN  FOR  GSS96  FOLLOWS  ',//! 

I F ( I COUNT . EQ« l IGQ  TO  385 

STOP 

END 
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